Application of the geometric measure of entanglement to three-qubit mixed states 
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The geometric measure of entanglement, originated by Shimony and by Barnum and Linden, is 
determined for a family of tripartite mixed states consisting of aribitrary mixtures of GHZ, W, and 
inverted-W states. For this family of states, other measures of entanglement, such as entanglement 
of formation and relative entropy of entanglement, have not been determined. The results for 
the geometric measure of entanglement are compared with those for negativity, which are also 
determined. The results for the geometric measure of entanglement and the negativity provide 
examples of the determination of entanglement for nontrivial mixed multipartite states. 
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Introduction: In a recent Paper [lj, we have explored 
the issue of quantifying entanglement by invoking certain 
simple elements of Hilbert space geometry. In doing this, 
we have been following a path originated by Shimony Q 
and by Barnum and Linden Q. Evidence for the rea- 
sonableness of this measure has been obtained via the 
investigation of certain bipartite mixed states for which 
other measures of entanglement (such as entanglement of 
formation and negativity) can also be computed. 

One of the virtues of the geometric approach to en- 
tanglement is its straightforward adaptability to any ar- 
bitrary multipartite systems (of finite dimension). In 
Ref. [l|, we described a procedure for applying the geo- 
metric measure of entanglement (GME) to certain highly 
symmetric multipartites states, and illustrated this pro- 
cedure via various examples. The aim of the present 
Paper is to present the analytical computation of the 
geometric measure of entanglement for a family of tri- 
partite mixed states consisting of aribitrary mixtures of 
GHZ, W, and inverted-W states (these pure states be- 
ing defined below). For this family, other measures of 
entanglement — such as entanglement of formation and 
relative entropy of entanglement — have not been com- 
puted analytically. (An exception is the measure of en- 
tanglement known as negativity, which, as we shall see, 
is readily computable.) 

We are motivated to study the quantification of entan- 
glement in multipartite mixed states for the basic reason 
that entanglement has been identified as a resource cen- 
tral to much of quantum information processing Q. To 
date, progress in the quantification of entanglement for 
mixed states has resided primarily in the domain of bi- 
partite systems [£j . For multipartite systems in pure and 
mixed states the quantification of entanglement presents 
even greater challenges. 

Basic geometric ideas; formulation of GME: We be- 
gin by briefly reviewing the formulation of this geo- 
metric measure in both pure-state and mixed-state set- 
tings. (For a discussion of the physical meanings of 
this measure, see Ref. 6|.) Let us start with a multi- 
partite system comprising n parts, each of which can 
have a distinct Hilbert space. Consider a general, n- 



partite, pure state (expanded in the local bases {le^}): 

~ Spi---p„ XpiP2---Pn \ e Pl e P2 ' 

the closest separable pure state 



) . As shown in , 
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satisfies the stationarity conditions 
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in which the eigenvalue A S [—1,1] is associated with 
the Lagrange multiplier enforcing the constraint {<f>\4>) = 
1, and the symbol denotes exclusion. Moreover, the 
spectrum A's can be interpreted as the cosine of the angle 
between and \<j>); the largest, A max , which we call 
the entanglement eigenvalue, corresponds to the closest 
separable state. We shall adopt E sin 2 = 1 — A^ iax as our 
entanglement measure for any pure state |i/>); we shall, 
in the following, drop the subscript 'max'. 

Given the definition of entanglement for pure states 
just formulated, the extension to mixed states p can be 
built upon pure states and is made via the use of the 
convex hull construction (indicated by "co" ) , as was done 
for the entanglement of formation (see, e.g., Ref. 0). 
The essence is a minimization over all decompositions 
P = J2iPi IV>i)(V>j| mto P ure states: 



E{ P ) 



(coE sin 2)(p) = min V ftB^d^)). (3) 
{Pi^i}^ 1 



This convex hull construction ensures that the measure 
gives zero for separable states; however, in general it also 
complicates the task of determining mixed state entan- 
glement. As mentioned in the Introduction, the principal 
aim of the present Paper is to calculate the GME for a 
specific, nontrivial class of tripartite mixed states. 

Arbitrary mixture of GHZ, W, and inverted-W states: 
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FIG. 1: Entanglement vs. the composition of the pure state 
IV'O.y)). 




0.4 0.6 

y 

FIG. 2: Entanglement of the pure state \ip{x = 0,y)) = 



y/y\W) + VI — V |W) vs. y. This shows the entanglement 
along the y axis. 

Of these, the least entangled state for given [x, y) has all 
coefficients non- negative (up to a global phase), i.e., 



Our goal is to calculate the GME for states of the form 

p(x,y) = x\GHZ)(GHZ\+y\W)(W\ + (l-x-y)\W)(W\, 

(4) 

where x,y > and x + y < 1. The three relevant pure 
states arc defined via 



\GHZ) = (|000) + | 111 
\W) = (|001) + |010) 
\W) = (|110) + |101> 



(5a) 

|100))/V3, (5b) 
|011))/V3, (5c) 



and have the following basic features. For \GHZ), 1 000) 
and | 111) are the closest separable states, and for it 
E sin 2(GHZ) = 1/2. For the W [or inverted- W] state, 

(V^73|0) + ^T/3\l))® 3 [or (V^73|l> + is 
one of the closest separable states, and E s [ n 2(W) = 

Esta>(W) = 5/9 > E sin2 (GHZ). 

A further property of the mixed state p(x,y), which 
facilitates the computation of its entanglement, is a cer- 
tain invariance, which we now describe. Consider the 
local unitary transformation a single qubit: 



10} - |0), 
/|1>, 



1) 



(6a) 
(6b) 



where g = exp(27ri/3), i.e., a relative phase shift. This 
transformation, when applied simultaneously to all three 
qubits, is denoted by U}.. It is straightforward to see that 
p{x, y) is invariant under the mapping 



1 3 

P:p^-*£u kP Ul 



(7) 



k=l 



Thus, we can apply Vollbrecht- Werner technique to 
the compution of the entanglement of p(x, y). 

The first step of this task is to find the general form of 
the pure states that, under P, are projected to p(x,y). 
This is readily seen to be 



\^(x, y)) = V^\GHZ) + y/y\W) + y/l-x- y\W). (9) 

The entanglement eigenvalue of \ip(x,y)) can then be 
readily calculated (see Ref. Q for the strategy), and one 
obtains 



A (x,y) = ^ t2 ^ [^(l+t 3 )+^t+^/3(l~x-y)t 2 y 

(10) 

where t is the (unique) non-negative real root of the fol- 
lowing third-order polynomial equation: 



+ ^3{l-x-y){-t 3 + 2t) = 0. (11) 

Hence, the entanglement function for \ip(x,y)), i.e., 
E^{x,y) = 1 — A(x,y) 2 , is determined (up to root- 
finding) . 




fx e^ 1 1 GHZ) + yfy~ e^ 2 1 W) + y/l-x-y e i<fc \W). (8) 



FIG. 3: Entanglement of the pure state \ij)(x,y = 1 — x)) — 
y/x\GHZ) + \/l — x \W) vs. x. This shows the entanglement 
along the diagonal boundary x + y = 1. Note the absence of 
convexity near x = 1; this region is repeated in the inset. 

Recall that our aim is to determine the entanglement 
of the mixed state p(x, y). As we now know the entangle- 
ment of the corresponding pure state \ij)(x,y)), we may 
accomplish our aim by invoking a result due to Vollbrecht 
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FIG. 4: Entanglement of the pure state \ip(x, (1 — x)rj) = 

^|Gf/Z) + V / (l - x)r \W) + ^{1 - x)(l - r)\W) vs. x and r. 
Note the symmetry of the surface with respect with r — 1/2. 




FIG. 6: Entanglement of the pure states \ip(x,(l — x)r)) = 

^x~\GHZ) + ^{l-x)r\W) + ^ {I - x){l - r)\W) vs. x for 
various values of r (from the top: 0, 0.1, 0.2, 0.3, 0.5). This 
reveals the nonconvexity in x in the (approximate) interval 
[0.85,1]. 




FIG. 5: Entanglement of the pure states |-0 (as, (1 — x)rY) = 

^\GHZ) + ^{l-x)r\W) + ^(l-x)(l-r)\W) vs. r for 
various values of x (from the bottom: 0.8, 0.85, 0.9, 0.92, 
0.94, 0.96, 0.98, 1). This reveals the nonconvexity in r for 
intermediate values of x. 




FIG. 7: Entanglement of the mixed state p(x,y). 



and Werner [Sj, which immediately gives the entangle- 
ment of p(x, y) in terms of that of \ip(x, y)), via the con- 
vex hull construction: E p (x,y) — (coE^)(x,y). Said in 
words, the entanglement surface z = E p (x, y) is the con- 
vex surface constructed from the surface z = E^(x,y). 

The idea underlying the use of the convex hull this. 
Due to its linearity in x and y, the state p of Eq. (J3J can 
[except when (x, y) lies on the boundary] be decomposed 
into two parts: p(x,y) — pp(xi, j/i) + (l -p)p(x 2 , y 2 ) with 
the weight p and end-points (xi,yi) and (x2,2/a) related 
by p xi + (1 — p)x2 = x and p y\ + (1 — p)yi — y ■ Now, if 
it should happen that pE^{x\,y\) + (1 — p)E^{x^y%) < 
E^,(x,y) then the entanglement averaged over the end- 
points gives a value lower than that at the interior point 
(x,y); this conforms with the convex-hull construction. 

It should be pointed out that the convex hull should 
be taken with respect to parameters on which the density 
matrix depends linearly, as x and y do in the example 
of p(x,y). Furthermore, in order to obtain the convex 
hull of a function, one needs to know the global struc- 
ture of the function; in the present case, E^,(x,y). We 
note that numerical algorithms have been developed for 
constructing convex hulls [9| . 



As we have discussed, our route to establishing the en- 
tanglement of p(x, y) invloves the analysis of the entan- 
glement of \ip(x, y)), which we show in Fig.^ Although 
it is not obvious, the corresponding surface fails to be 
convex near to the point (x, y) = (1, 0), and therefore, in 
this region, we must suitably convexify in order to obtain 
the entanglement of p(x,y). To illustrate the properties 
of the entanglement of \ip(x,y)) we show, in Fig. the 
entanglement of \ip(x, y)) along the line x = 0; evidently 
this is convex. By contrast, along the line x+y = 1, there 
is a region in which the entanglement is not convex, as 
Fig. 13 shows. The nonconvexity of the entanglement of 
\tp{x, y)) complicates the calculation of the entanglement 
of p(x,y), as it necessitates a procedure for constructing 
the convex hull in the (as it happens, small) nonconvex 
region. Elsewhere in the xy plane the entanglement of 
p{x, y) is given directly by the entanglement of \tp(x, y)). 

At worst, convexification would have to be undertaken 
numerically. However, in the present setting it turns out 
one can determine the convex surface essentially analyt- 
ically, by performing the necessary surgery on surface 
z = E^(x, y). To do this, we make use of the fact that if 
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we parametrize y via (1 — x)r, i.e., we consider 

p(x,(l-x)r) = x\GHZ)(GHZ\ + (l-x)r\W)(W\ 
+ (l-x)(l-r)\W)(W\, (12) 

where < r < 1 [and similarly for \ip(x,y)}] then, as 
a function of (x,r), the entanglement will be symmetric 
with respect to r = 1/2, as Fig. 0] makes evident. With 
this parametrization, the nonconvex region of the entan- 
glement of \ip) can more clearly be identified. 

To convexify this surface we adopt the following conve- 
nient strategy. First, we reparametrize the coordinates, 
exchanging y by (1 — x)r. Now, owing to the linearity, in 
r at fixed x and vice versa, of the coefficients x, (1 — x)r 
and (1 — x) (1 — r) in Eq. I|12fl . it is certainly necessary for 
the entanglement of p to be a convex function of r at fixed 
x and vice versa. Convexity is, however, not necessary in 
other directions in the xr plane, owing to the nonlincar- 
ity of the the coefficients under simultaneous variations 
of x and r. Put more simply: convexity is not neces- 
sary throughout the (x, r) plane because straight lines in 
the (x, r) plane do not correspond to straight lines in the 
(x, y) plane (except along lines parallel either to the r or 
the x axis). Thus, our strategy will be to convexify in 
a restricted sense: first along lines parallel to the r axis 
and then along lines parallel to the x axis. Having done 
this, we shall check to see that no further convexification 
is necessary. 

For each x, we convexify the curve z = (x, (1 — x)r) 
as a function of r, and then generate a new surface by 
allowing x to vary. More specifically, the nonconvexity in 
this direction has the form of a symmetric pair of minima 
located on either side of a cusp, as shown in Fig.0 Thus, 
to correct for it, we simply locate the minima and connect 
them by a straight line. 

What remains is to consider the issue of convexity 
along the x (i.e., at fixed r) direction for the surface just 
constructed. In this direction, nonconvexity occurs when 
x is, roughly speaking, greater than 0.8, as Fig. sug- 
gests. In contrast with the case of nonconvexity at fixed 
r, this nonconvexity is due to an inflection point at which 
the second derivative vanishes. To correct for it, we lo- 
cate the point x = xq such that the tangent at x = xq is 
equal to that of the line between the point on the curve at 
Xq and the end-point at x — 1, and connect them with a 
straight line. This furnishes us with a surface convexified 
with respect to x (at fixed r) and vice versa. 

Armed with this surface, we return to the (cc, y) 
parametrization, and ask whether or not it is fully con- 
vex (i.e., convex along straight lines connecting any pair 
of points). Said equivalently, we ask whether or not any 
further convexification is required. Although we have 
not proven it, on the basis of extensive numerical ex- 
ploration we are confident that the resulting surface is, 
indeed, convex. The resulting convex entanglement sur- 
face for p(x, y) is shown in Fig. [7| Figure [S] exemplifies 
this convexity along the line 2y + x = 1. We have ob- 
served that for the case at hand it is adequate to correct 
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FIG. 8: Entanglement of the mixed state 

p(x,V=(l-x)/2) = x\GHZ)(GHZ\ + 1=*(\ W )(W\ + 
\W){W\) vs. x. Inset: enlargement of the region x £ [0.2,0.3]. 
This contains the only point, (x,y) = (1/4,3/8), at which 
E p (x,y) vanishes. 

for nonconvexity only in the x direction at fixed r. 

Negativity: This measure of entanglement is defined as 
twice the absolute value of the sum of the negative 
eigenvalues of the partial transpose of the density ma- 
trix 0, ^J. In the present setting, viz., the family 
p(x, y) of three-qubit states, the partial transpose may 
equivalently be taken with respect to any one of the three 
parties, owing to the invariance of p(x, y) under all per- 
mutations of the parties. Transposing with respect to the 
third party, one has 

N(p) = ~2J2 K (13) 

Ai<0 

where the A's are the eigenvalues of the matrix p Tz , 

It is straightforward to calculate the negativity for 
p(x,y); the results are shown in Fig. It is interest- 
ing to note that, for all allowable ranges of {x,y), the 
state p(x,y) has nonzero negativity, except at {x,y) — 
(1/4,3/8), at which the calculation of the GME shows 
that the density matrix is indeed separable. The fact 
that the only positve-partial-transpose (PPT) state is 
separable is the statement that there are no entangled 
PPT states (i.e., no PPT bound entangled states) within 
this family of three-qubit mixed states. The negativity 
surface, Fig. [51 is qualitatively — but not quantitatively — 
the same as that of GME. By inspecting the negativity 
and GME surfaces one can see that they present ordering 
difficulties. We mean by this that one can find pairs of 
states p(x\,yi) and p(x2,y2) that respectively have neg- 
ativities Ni and N% and GMEs E\ and E2 such that, say, 
Nt < N2 but Ei > E2- Said equivalently, the negativity 
and the GME do not necessarily agree on which of a pair 
of states is the more entangled. For two qubit settings, 
such ordering difficulties do not show up for pure states 
but can for mixed states [TTl[T^ | . On the other hand, such 
difficulties already show up for pure states, as the follow- 
ing example shows: N(GHZ) = 1 > N(W) = 2\/2/3 
whereas for the GME the order is reversed. We note 
that for the relative entropy of entanglement Ejj, one 
has E R (GHZ) = log 2 < E R (W) = log(9/4)[l3|. 
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FIG. 9: Negativity of the mixed state p(x,y). 



Concluding remarks: By making use of the geometric 
measure of entanglement we have addressed the entan- 
glement of a rather general family of three-qubit mixed 
states analytically (up to root-finding). This family con- 
sists of arbitrary mixtures of GHZ, W, and inverted-W 
states. To the best of our knowledge, corresponding re- 
sults have not, to date, been obtained for other measures 
of entanglement, such as entanglement of formation and 



relative entropy of entanglement. We have obtained cor- 
responding results for the negativity measure of entan- 
glement, and have compared them with those for the ge- 
ometric measure of entanglement. Among other things, 
we have found that there are no PPT bound entangled 
states within this general family. 

We are unaware of any explicit generalization of entan- 
glement of formation to multipartite mixed states. How- 
ever, if such a generalization should emerge, and if it 
should be based on the convex hull construction (as it 
is in the bipartite case), then one may be able to cal- 
culate the entanglement of formation for the family of 
mixed states considered in the present Paper. It would, 
then, be interesting to know whether or not the similari- 
ties between entanglement of formation and the geomet- 
ric measure of entanglement found at the level of certain 
bipartite mixed states y| continue to hold beyond the 
bipartite world. 
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